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SUMMARY 

Evaluation is made of two computer programs capable of predicting the 
potential and viscous interacting flow around wings of infinite aspect ratio. 

The programs are compared in terms of their capabilities, the approximations 
and the methods of solution used, and the input requirements. Six airfoils, 
each representative of a class of airfoils, are used as test airfoils. The 
results predicted by the programs are presented for each airfoil at sweep angles 
of 0°, 20®, and 40® over a range of angles of attack. The results show that at 
zero sweep both programs predicted the aerodynamic coefficients well and 
generally in good agreement with measurements. At 20® and 40® of sweep, as 
there are no experimental data available, definitive conclusions cannot be 
drawn about the accuracy of the predictions although the results are presented 
and discussed. The execution times are approximately the same for the two 
programs. 


INTRODUCTION 

The present day aerodynamics designer has a number of computer programs 
(ref. 1) available as tools for design. While most of these programs are based 
on potential -flow analysis, a limited number (refs. 2 and 3) are available that 
solve the potential -viscous interacting flow in two dimensions. However, to 
date the author knows of only two programs that predict the aerodynamic charac- 
teristics of swept- infinite wings. (See refs. 4 and 5.) Swept-infinite wings 
are swept wings of infinite aspect ratio or infinite span. Although swept- 
infinite wings do not appear to be directly usable in aerodynamics or hydro- 
dynamics, knowledge of * the prediction methods for swept-infinite wings should 
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be useful and perhaps the necessary first step in the ultimate goal of the 
prediction of the characteristics of finite wings. 

The aim of this paper is to compare, from the user's viewpoint, the swept- 
infinite wing potential/viscous flow programs of Dvorak and Woodward (ref. 4) 
anl Gingrich and Bonner (ref. 5) for single-element wings. Because of the 
proprietary nature of the latter program, details of it will be discussed only 
sketchily here. However, six test airfoils are run on both programs, and the 
results are presented and compared. In the unswept orientation, the two predic- 
tions are also compared with those of reference 3 with some in-house modifi- 
cations made at NASA Langley Research Center (ref. 6). 

SYMBOLS 

The units used for the physical quantities of this paper are given both in 
the International System of Units (SI) and in the U.S. Customary Units. 

The symbols enclosed in parenthesis are used in plotting the results. 


c 

airfoil normal chord, cm (in 

.) 

Cd(CD) 

drag coefficient. 

D 

1 II 


Cl(CL) 

lift coefficient. 

L 

1 II 2r 

7 ^ 


D 

drag force per unit span, N 

(Ibf) 

L 

lift force per unit span, N 

(Ibf) 

M 

on 

free-stream Mach 

number 


RN 

Reynolds number = 

U C/v 

oo ' 


u. 

free-stream velocity, m/sec 

(ft/sec) 


2 


X coordinate along the normal chord line, cm (ft) 

a(ALPHA) angle of attack (angle between the streamwise chord and the free- 
stream direction), deg 

A sweep angle, deg 

2 2 

V ‘ kinematic viscosity of air, m /sec (ft /sec) 

Q 3 

p density of air, kg/m (slugs/ft ) 

COMPARISON OF PROGRAMS 


The two programs are compared In terms of their capabilities, the approxi- 
mations used, and the methods of solution In the following table: 


DVORAK AND WOODWARD PROGRAM (DWP] 


GINGRICH AND BONNER PROGRAM 


1. Program capabilities: 

(a) Has two-dimensional and - Same - 

swept capability. 

(b) Has multi -element capability Has capability for a single element 

up to one slat, one main sharp trai ling-edge airfoil only, 

element and two flaps. 

2. Overall method of analyses: 

Analysis Is based on potential/ - Same - 
viscous Interacting flow. 

3. Type of flow assumed: 

Incompressible flow analysis Incompressible flow analysis with 

Is used. Labrujere (ref. 7) corrections for 

compressibility effects Is used. 

4. Simplifications used in the 3-D flow: 

Independence principle Is used. - Same - 

5. Method of potential -flow solution used: 

Pressure distribution Is deter- Pressure distribution Is determined 
mined perpendicular to the leading using the two-dimensional conformal 
edge using a two-jdimenslonal vor- transformation of Theodorsen for 
tex sheet technique with an Incompressible flow, 

alternate formulation for blunt 
trailing edges. 
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OWP 


GBP 


6. Coupling between potential and viscous solutions: 


Boundary- layer displacement 
effect Is Introduced Into the 
potential flow by an equivalent 
source distribution. Thus elimi- 
nating the necessity of Inverting 
a new Influence coefficient matrix 
for each inviscld-viscous 
iteration. 

7. Type of boundary- layer method used: 

Integral boundary- layer methods 
are used on all surfaces except 
the upper surfaces of the flaps. 


The method of Spence (ref. 8) Is 
used to establish the effect of the 
displacement surface on the potential 
flow. 


Integral boundary- layer method are 
used on both surfaces of the single 
element. 


8. Determination of program starting conditions: 

To start computations, integral - Same - 
properties on the stagnation line 
are determined by the method of 
Cumpsty and Head (ref. 9). 


9. Method of laminar boundary- layer solution: 


Laminar boundary- layer growth is 
calculated using a two-dimensional 
Integral method along external 
streamlines. A check is made for 
laminar separation, free transi- 
tion, and forced transition. 

10. Determination of laminar separation: 


11. Determination of transition: 

Transition criteria are based on 
the streamwise boundary- layer 
parameters (refs. 11 to 13). 
Crossflow Instability is not 
considered. 


For laminar boundary- layer calcula- 
tions, the approximate three- 
dimensional Integral method of Smith 
and Young (ref. 10) Is used with 
extensions to adiabatic compressible 
flow. 


Separation Is predicted when shear 
stress in the normal chord direction 
becomes zero. 


- Basically the same (see refs. 11 
and 14) - 

Crossflow instability analysis of 
Brown (ref. 15) is extended and used. 


A correlation between the pressure 
gradient parameter and Reynolds 
number based on momentum thick- 
ness Is used to predict laminar 
separation and possible reattach- 
ment as a turbulent boundary layer. 
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12. Method of turbulent boundary-layer solution: 

Turbulent boundary- layer develop- Turbulent boundary layer Is treated 
ment In the streanwlse and the with the three-dimensional Integral 
crossflow direction Is calculated method of Smith (ref. 18). Small 

by the three-dimensional Integral crossflow assumption Is not made, 

method of Cumpsty and Head (ref. 

16). The relationship between the 
streamwise and crossflow profiles 
suggested by Mager (ref. 17) Is 
also used. 

13. Determination of turbulent separation: 

Turbulent separation Is predicted Separation is predicted when shear 

when the angle between the normal stress In the normal chord direction 

chord and skin-friction direction becomes zero, 

becomes 90*. 

14. Determination of drag coefficient: 

Drag coefficient Is determined Squire and Young formula Is extended 
from Squire and Young (ref. 19) for small crossflow and used to 
formula using streamwise boundary- determine the drag coefficient, 
layer parameters at the trailing 
edge. 

15. Convergence of the potential -viscous solution: 

During any Iteration, if the lift Convergence criterion is that the 

coefficient Is within 0.01-0.015 change In circulation between two 

(the actual in-between value to successive iterations be less than a 

be chosen by the user) of the certain value. This condition Is 

previous Iteration, the solution enforced when the difference between 
Is considered to converge. the velocities on the upper and lower 

surfaces become a small fraction of 
the free-stream velocity. 


INPUT-OUTPUT REQUIREMENTS 


The purpose of this sec^'ion is to compare from the user's viewpoint the two 
programs In terms of their input and storage requirements. For convenience, the 
comparisons are tabulated: 
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DV ORAK AND WOODWARD PROGRAM (DWP) 

1. Airfoil geometry is to be Input by 
the coordinates of upper and lower 
surfaces. In each case, from the 
L.E. to the T.E, .Intersection of 
the surface with the reference line 
i.xjst be the L.E. point (first point 
on each surface). The coordinates 
should be Input dimensionless but 
not In percentages. 

2. Maximum number of Input coordinate 
points is limited to 30 on each of 
the surfaces. 

3. Flow parameters to be input are the 
free-stream RN, Mach number and 
velocity, and the normal chord 
length. For single-element cases, 
the program uses only the RN. 

4 . Program interprets the input angle of 
attack as the angle relative to the 
Input reference chord line. 


GINGRICH AND BONNER PROGRAM (GBP) 

Coordinates of the wing are to be 
Irput from the leading edge (L.E.) 
to the trailing edge (T.E.) on the 
upper surface and continued from the 
T.E. to the L.E. on the lower surface 
In a clockwise direction. 


Total number of Input points around 
the airfoil is limited to 50. 


Flow parameters to be input are the 
free-stream Mach number, stagnation 
temperature and pressure and the 
normal (not streamwise) chord. 

Program calculates the RN. 

Program calculates the longest chord- 
line and the Input angle of attack Is 
taken with respect to the longest 
chord line. 


NOTE; In both programs, the input parameters should be based on the 
geometric parameters in the unswept orientation of the \/1ng. 


5. Transition can be free or forced. 
For forced transition, x/c value at 
which transition Is desired should 
be Input. Both the upper and lower 
surfaces get tripped at the same 
x/c location. 

6. In the event of separation without 
reattachment, the output does not 
clearly state the result. Instead 
the computations continue without 
interruption. It appears that the 
condition for separation without 
reattachment is not satisfied for 
any of the six test airfoils. 

7. Aerodynamic coefficients are based 
on the chord length normal to the 
L.E. 


Transition can be free or forced. 

For forced transition, the input 
coordinate point number that repre- 
sents the x/c at which transition Is 
desired should be specified. 


The program does not model the possi 
blllty of reattachment in the event 
of laminar separation. 


Aerodynamic coefficients are based on 
the streamwise chord length. 
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DWP 


GBP 


8. Program required 150K octal words 
of memory. 

9. A numerical value for convergence 
criterion on C|_ Is required as 
an Input. The user can choose a 
value between 0.01 and 0.015. The 
user can also Input the maximum 
number of Iterations as an addi- 
tional parameter. 


It required under 65K octal words 
of memory . 

User has no direct control on 
convergence criterion. Indirectly, 
the user may Input the maximum 
number of iterations as a parameter. 
This number is not to exceed 15. 


RESULTS AND DISCUSSION 


'he following six airfoils, each representative of a class of airfoils, as 
Indicated, were chosen as the test cases: (1) NACA 0012 airfoil, uncambered, 

12 percent ^hlck; (2) NACA 2424 airfoil, slightly cambered but with a much 
higher thickness of 24 percent chord; (3) NACA 63-006 airfoil, symmetric 6- 
percent-thlck airfoil designed for low drag; (4) NACA 23012 airfoil, cambered 
version of NACA 0012 airfoil; (5) 6A(W)-1 airfoil, advanced technology general 
aviation airfoil designed for wide C|_ range without flow separation, 17- 
percent-thlck with substantial camber near the trailing edge, a reflexed lower 
surface contour, and a blunt trailing edge; and (6) TN-D 7071 airfoil (see ref. 
24), designed for maximum lift, maximum thickness 12.5-percent chord, located at 
approximately 25 percent chord from the leading edge. 

For each airfoil, one computer job was run at a fixed sweep angle In order 
to compute a number of angle-of-attack cases. Three such jobs, one each corre- 
sponding to the sweep angles of 0®, 20®, and 40®, were run for each airfoil 
through the two programs under consideration. For the unswept orientation, the 
airfoils were also run through the program of reference 3 updated with In-house 
modifications at Langley (ref. 6). 

In the unswept orientation for each of the test-case airfoils, aerodynamic 
coefficients predicted by the three programs are compared with the experimental 
results available In the literature. For the Infinite-span airfoils In the 
swept orientation, experimental data for C|^ and Cq do not exist In the open 
literature at the present time. There Is some published boundary- layer data for 
a few swept-wing cases (refs. 20 and 24) primarily for the purpose of evaluating 



three-dimensional turbulent boundary- layer computational methods. Therefore, 
conclusions cannot be drawn about agreement between the predictions and 
measurements for wings In swept orientation. 

The airfoil -surface coordinates, chord length, Reynolds number, and the 
angle of attack Input Into both programs are based on the unswept orientation 
or what Is called by the authors of both programs as the reference orientation. 
However, the output aerodynamic coefficients of Dvorak and Woodward are based 
on the reference (that Is, normal) chord; whereas, those of Gingrich and Bonner 
are based on the streamwise chord. It Is, therefore, necessary to point out 
that all the results presented In this report are based on the streamwise chord. 

There Is one other point that requires mention here. For all cases, the Gingrich 
and Bonner Program (GBP) first calculates the longest chord and the input angle 
of attack Is then Interpreted to be the angle of attack of the flow with respect 
to the longest chord. For most airfoils, the reference line (the line with 
respect to which the airfoil surface coordinates are specified) is coincident 
with the longest chord line. Of the six test-case airfoils, the TN D-7071 
airfoil Is the only one with a reference line that Is not coincident with the 
longest chord line. The angle between the two lines for this airfoil is 5.03476®. 
However, experimental as well as the Dvorak and Woodward results for all airfoils 
are based on the angle of attack defined with respect to the Input reference 
line. Therefore, for proper comparisons with these results, the angles of 
attack for the TN D-7071 airfoil are adjusted by adding 5.03476° to the value 
of the angle of attack listed on the output. 

Both programs have the ability to continue computations, without encountering 
numerical difficulties, for as much as 15° (depending on the airfoil) past the 
angle of attack at which trailing edge separation has first occurred. In other 
words, calculatations can be continued beyond the case of the small separation 
by incrementing the angle of attack until the boundary layer has separated over 
much of the upper surface even though such predictions cannot be meaningful 
because the boundary- layer approximations used in the analysis are no longer 
valid. Therefore, results are presented only for those angle-of-attack cases In 
which there Is no separation at all or separation is close to the tral ling-edge 
region within 10 percent of the chord length from the trailing edge. 
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For the NACA 0012 airfoil at zero sweep, all three programs predict 
and Cq values that are In good agreement with experimental data as shown In 
figure 1(a). At higher angles of attack, the Dvorak and Woodward Program (DWP) 
underpredicted C|_ while the (>BP overpredicted C|_. When these C|^ predictions 
trends are taken Into account, DWP and GBP underestimate the drag coefficient at 
all angles of attack »whereas the Modified NASA Program (MNP) overestimated It. 

At 20^ sweep, both DWP and GBP predict the same values for C|^, but relative to 
GBP, DWP predicts higher drag coefficients. 

For the NACA 2424 airfoil In the unswept orientation, figure 2 shows that 
all three programs predicted linear C|_ versus a curves. Ci_ values of DWP 
are In s^greement with MNP and both of these predictions were consistently higher 
than thDse of GBP. The experimental Cl versus a curve is nonlinear pri- 
marily due to the higher thickness of the airfoil, but all these programs falle . 
to predict this nonlinear thickness effect. In view of the Cl predictions, 

MNP overpredicted Cq; where=*s DWP and GBP underestimated Cq. At 20® and 40° 
of sweep, DWP consistently predicted higher Cl and Cq values over GBP. 

For the NACA 63-006 airfoil, figure 3 shows that for the unswept orientation, 
predictions of all programs are generally In good agreement with the experimental 
results. In the swept orientation, predictions of both DWP and GBP are In close 
agreement. 

For the NACA 23012 airfoil, all three programs predicted (see fig. 4) linear 
Cl versus a curves for the unswept wing. DWP and MNP predictions of Cj^ are i 
good agreement with experimental data, while GBP predicted consistently higher 
values for Cl. The predictions of Cq by MNP are in good agreement with 
measurements. Figure 4(b) shows that DWP and GBP have underpredicted Cq, and 
the extent of disagreement with measurements can be seen to increase with 
Increasing angle of attack. For the 20° swept case. Cl predictions by GBP are 
higher than those by DWP; whereas, DWP predicted higher values for Cq than did 
GBP. Further DWP predicted smaller stall angle. At 40° sweep, among the two 
programs DWP predicted higher values for both Cl and Cq. 

The results for the TN D-7071 airfoil are shown in figure 5. The C^^ and 
Cq predictions by all the three programs are In better agreement for this air- 
foil than for any of the other airfoils investigated in this report. Such good 





9 



agreement Is rather surprising considering the fact that this Is an unusual air- 
foil designed for high lift. Departures In predictions occurred very close to the 
stall region. No conclusions can be drawn on Cq predictions because there are 
no experimental data available on Cq even for the 0** sweep case. At 20° sweep, 
the results are much like those at 0* sweep. At 40° sweep. GBP predicted early 
separation while DUP did not. 

The results for the GA(W)-1 airfoil are shown In figure 6. At 0° sweep, 
predictions of by DWP and MNP are In good agreement with experiment; where- 
as. GBP consistently underpredicted C|^. Although both DWP and GBP underpredicted 
Cq • the predictions are In good agreement with each other. MNP slightly under- 
predicted Cq. To eliminate the possible differences In results due to transi- 
tion modeling, this airfoil was run at 0° sweep with fixed transition. The 
results are shown In figures 6(a) and 6(b). For each program, the trends 
predicted with free and fixed transition a'^e similar and do not appear to be 
effected significantly by the transition modeling. At 20® and 30° sweeps. DWP 
predicted consistently higher values for C|_ than did GBP. Both of these 
programs predict about the same values for Cq. 

The DWP was run In all cases with the parameters for maximum number of 
Iterations set at 6 and the value for convergence criterion on C|_ set at 0.01. 
For the six airfoils tested, the convergence criterion was observed to be satis- 
fied In under six Iterations In most of the cases and only in a few cases th» 
sixth Iteration was caaputed. In all the rases where the sixth Iteration was 
computed. C|_ was observed to be wUh<n 0.015 of the previous Iteration. The 
above comments are applicable only In the absence of catastrophic separation. 

The GBP was run In all cases with the parameter for maximum number of 
iterations set at seven. C^^ Is observed to converge In all cases to a value 
within 0.01 of the previous Iteration In the absence of catastrophic separation. 

CONCLUSIONS 

Six single-element airfoils are used as test cases to compare the perfor- 
mance of the throe programs by Dvorak and Woodward (ref. 4). Gingrich and Bonner 
(ref. 5). and the modified NASA program (refs. 3 and 6). Based on the results 
obtained, the following conclusions can be drawn: 
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1. In the unswept orientation of the airfoils, all the three programs 
predicted the lift coefficient well for all the "rfolls except for the 
percent- thick airfoil. Further, the predicted versus a curves are approxi- 
mately linear for all the airfoils. In the swept orientation, the predictions 

of C|^ by tno two swept-wing programs are In good agreement for the NACA 0012, 

NACA 63-006, and I'N 0-7071 airfoils; whereas, for the NACA 2424 and GA(U)-1 air- 
foils, the Dvorak and Woodward program consistently predicted higher values 
for C^. 

2. For the airfoil designed for high lift, somewhat surprisingly all the 
programs predicted Cj_ accurately, and while there are no experimental data 
available to compare the drag coefficients with, the Cp predictions of all 
programs are In good agreement with each other. 

3. In the unswept orientation, the drag coefficients are predicted satis- 
factorily at small angles of attack, and the predictions by the swept-wing programs 
progressively underestimated the drag coefficient as tne stall region Is approached. 
In general, the predictions of both Cl and Cp by the modified NASA program 
(which has no sweep capability) are in better agreement with measurements than 

the two swept-wing programs. For the airfoils in the swept orientation, as 
there are no published measurements, definitive conclusions cannot be drawn 
about the Cp predictions. At 20° sweep, Dvorak and Woodward program tended 
to predict higher Cp values than the Gingrich and Bonner program. However, at 
40° sweep, there Is little agreement between the predictions of the two programs. 

4. The execution times are approximately the same for all the three 
programs. In general, the time increased with Increasing sweep angles. For an 
angle-of-attack case, the typical execution times on the COC 6600 computer are 
20, 25, and 30 seconds, respectively, corresponding to the 0°, 20°, and 40° sweep 
angles. For a given case, the execution times of the programs are typically 
within 5 seconds of each other. The memory required for the Dvorak and Woodward 
program Is approximately 150K octal words and that for the Gingrich and Bonner 
program Is approximately 65K. The former is higher primarily because of the 
multi -element capability of that program. 
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igure 3. - Comparison of predicted lift and drag ocefficients of NACA 63^006 airfoil 
at RN ■ 3.335 x 10^, M ■ 0.14, and with free transition. 


0 Dvorak & Woodward 
^ Gingrich & Bonner 



Figure X - Continued. 
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Figure 4. - Continued. 
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Fiaure 4. - Concluded 
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Figure 5. - Comparison of predicted lift and drag coefficients of TN D-7071 airfoil 
at RN - 9.06 X #, M - 0.14, and with free transition. 
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(c) A* 20 , free transition 
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Figure 6. - Concluded 



